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ABSTRACT 

• »"H . Aims. We present results from deep y-ray observations of the Galactic pulsar wind nebula HESS J1825— 137 performed with the H.E.S.S. array. 
' Methods. Detailed morphological and spatially resolved spectral studies reveal the very high-energy (VHE) y-ray aspects of this object with 
, unprecedented precision. 

■ - - Results. We confirm previous results obtained in a survey of the Galactic Plane in 2004. The y-ray emission extends asymmetrically to the 
south and south-west of the energetic pulsar PSR J1826-1334, that is thought to power the pulsar wind nebula. The differential y-ray spectrum 
of the whole emission region is measured over more than two orders of magnitude, from 270 GeV to 35 TeV, and shows indications for a 
deviation from a pure power law. Spectra have also been determined for spatially separated regions of HESS J1825-137. The photon indices 
from a power-law fit in the different regions show a softening of the spectrum with increasing distance from the pulsar and therefore an energy 
dependent morphology. 

Conclusions. This is the first time that an energy dependent morphology has been detected in the VHE y-ray regime. The VHE y-ray emission 
of HESS J1825— 137 is phenomenologically discussed in the scenario where the y-rays are produced by VHE electrons via Inverse Compton 
scattering. The high y-ray luminosity of the source cannot be explained on the basis of constant spin-down power of the pulsar and requires 
higher injection power in past. 

Key words. ISM: plerions - ISM: individual objects: PSRB1823-13, HESS J 1825-1 37, G18.0-0.7 - gamma-rays: observations 



1. Introduction 

A growing number of extended objects that seem to be as- 
sociated with energetic pulsars are detected in the Galactic 
Plane by their very high-energy (VHE, energy E Y > 
100 GeV) y-ray emission. Latest results on this class of 
object s include emission from MSH-15-52 (HESS J1514- 
591) JAharonian et all l2005al) and Vela X (HESS J0835- 
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455) JAharonian et alJ l2006al) . and the two sources in the 
Kookaburra regi on (HESS J 1420-607 and HESS J 14 18-609) 
as described in lAharonian et al.l d2006bl) . If these associa- 
tions are correct, then these objects are pulsar wind nebulae 
(PWN), objects generally thought to be powered by a rela- 
tivistic particle outflow (electrons and positrons) from a cen- 
tral source. The central source - a pulsar - is a rapidly rotat- 
ing neutron star generated in a supernova event. The relativis- 
tic wind of particles flows freely out until its pressure is bal- 
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Fig. 1. Acceptance-corrected smoothed excess map (smoothing radius 2.5') of the 2.7° x 2.7° field of view surrounding 
HESS J1825— 137. The linear colour scale is in units of integrated excess counts within the smoothing radius of 2.5'. The ex- 
cess has been derived from a model of the system acceptance as described in the text. The inset in the bottom left corner shows 
the PSF of the dataset (smoothed in the same way as the excess map with the black dashed circle denoting the smoothing radius). 
The dashed black and white contours are linearly spaced and denote the 5cr, 10cr and 15cr significance levels (the 5<x contour 
being the outermost one), determined with a 9 — 0.1° radius cut. The best fit position of HESS J1825-137 is marked with a 
black square, the best extension and position angle by a black ellipse (see text). The dotted white contour shows the 95% po- 
sitional confidence contour of the unidentified EGRET source 3 EG J 1826-1 302. The position of the pulsar PSR J 1826-1 334 is 
mark ed by a white triangle. T he bright point-source to the south of HESS J 1825-1 37 is the microquasar LS 5039 (HESS J 1826- 
148) JAharonian et all2005dh . The colour scale for this source is truncated in this Figure. The Galactic plane is shown as a white 
dashed line. Some indication for an additional emission region to the north of the pulsar is seen. 



anced by that of the surrounding medium. In that region the 
wind decelerates and a standing te rmination shock is formed 
at which particles are accelerated jKennel & Coronitilll984al 
lAharonian. Atovan & Kifunel 19971) . The existence of electrons 
accelerated to energies > 100 TeV in such PWN has been es- 
tablished by X-ray observations of synchro tron emission, e.g. 
in the Crab nebula JWeisskopf et aimOOOh . VHE y-rays can 
be generated in PWN from the high-energy electrons by non- 
thermal bremsstrahlung or inverse Compton (IC) scattering on 
photon target fields, such as the cosmic microwave background 
(CMBR) or star-light photons. 

One such object, HESS J1825-137, has been detected 
by the High Energy Stereo scopic System (H.E.S.S.) in a 
survey of the inner Galaxy ( lAharonian et alJ l2005bl l2006ci) 
and has subsequently been associated with the X-ray PWN 
G18.0 -0.7 surrounding the energetic pulsar PSRJ1826- 
1334 (lAharonian et alJ l2005ch . This pulsar PSR J 1826-1 334 
(also known as PSR B 1 823-1 3) was detect ed in the Jodrell 
Bank 20 cm radio survey IClifton et alJ 19921) and is among the 



20 most energetic pulsars in the current ATNF catalogue (spin 
down power £ = 3 x 10 36 erg/s). The distance of PSR J1826- 
1334 as meas ured from the dispersi on of the radio pulses is 
3.9 + 0.4 kpc JCordes & Laziol2002h . The radio detection fur- 
ther revealed characteristic properties of the system that are 
similar to those of the well studied Vela pulsar, namely a pulse 
period of 101 ms and a characteristic age of 21.4 kyears (de- 
rived by t = P/2P). This age renders PS R J 1826- 13 34 one 
of the 40 youngest pulsars detected so far jManchester et alJ 
I2005I) . and due to this, deep radio observations were performed 
to find emission associated with the remnant of the Supernova 
explosion that gave rise to the pulsar. However, deep VLA 
observations of the 20' surrounding th e pulsar have failed 
to det ect this Supernova remnant (SNR) jBraun. Goss & Lvnel 
119891) . 

Initi al observations o f the region in X-rays with 
ROSAT frinlev et al J 19981) revealed a point source surrounded 
by an elongated diffuse region of size ~5'. The X-ray emis- 
sion region was subsequentially observed with the ASCA in- 
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strument and the data confirmed the picture of a compact ob- 
ject su rrounded by an extended emission region dSakurai et alJ 
l200ll) . While ROSAT data did not provide sufficient statistics, 
ASCA data lacked the spatial resolution to resolve and inter- 
pret the sources in this region. The situation was clarified in 
an XMM-Newton observation in which high angular resolu- 
tion observations revealed a compact core of extension 30" 
surrounding PSR J 1826-1 334, and furthermore an asymmetric 
diff use nebula extending at least 5' to the south of the pul- 
sar JGaensler et alJl2003l) . In this XMM-Newton dataset the 
signal to noise ratio deteriorates rapidly at offsets larger than 
5' and for this reason the XMM data cannot place useful con- 
straints on the presence of a faint shell of emission at larger 
radii as might be produced by an associated SNR. The extended 
asymmetric structure was attribute d to synchrotron emi ssion 
from the PWN of PSR J 1826-1 334 JGaensler et alJ2003l) . The 
X-ray spectrum in the diffuse emission region follows a power 
law with photon index F ~2.3 and an X-ray luminosity between 
0.5 and 10 keV of L x ~ 3 x 10 33 erg s -1 compared to the X- 
ray spectrum for the compact core following a power law with 
F ~1.6 and L x ~ 9x 10 32 erg s ' 1 (these luminosities a re derived 
assuming a distance of 4 kpcV lGaensler et all J2003h discussed 
various scenarios to explain the asymmetry and offset morphol- 
ogy of the PWN G18.0-0.7. The most likely explanation seems 
to be that a symmetric expansion of the PWN is prevented by 
dense material to the north of the pulsar which shifts the whole 
emission to the south. Asymmetric reverse shock interactions 
of this kind have originally been proposed to explain the off- 
set morphology of the Vela X PWN based on h y dro-dy namical 
simulations by iBlondin. Chevalier & Friersonl J200 lh . Indeed 
recent analyses of CO data show dense material surrounding 
PSR J1826-1334 (at a dist ance of 4 kpc) to the north and north- 
east JLemiere et al.ll20"05h . supporting this picture. It is inter- 
esting to note, that H.E.S.S. has now detected offset morpholo - 
gies from both G18. 0-0.7 and Vela X JAharonian et alJ2006al) . 
confirming the existence of a c lass of at least two o ffset PWN 
implied by X-ray observations JGaensler et alJ2003l) . Whereas 
X-rays probe a combination of the thermal and ultrarelativistic 
components, which could have been mixed at the time when the 
asymmetric reverse shock interaction took place, the H.E.S.S. 
results are important in determining the offset morphology of 
the ultrarelativistic component alone. 

Based on its proximity and energetics, the pulsar 
PSR J 1826-1 334 has been proposed to be asso ciated with the 
unide ntified EGRET source 3EG J1826-1302 JHartman et all 
Il999l) . This EGRET source exhibits a hard power law of pho- 
ton index 2.0 + 0.11 with no indication of a cut-off. The pul- 
sar lies south of the centre of gravity of the EGRET position 
and is marginally enclosed i n the 95% confidence c ontour (see 
Fig.Q. It has been shown JZhang & Chengiri998h that an as- 
sociation between PSR J1826-1334 and 3EG J1826-1302 is 
plausible based on the pulsar properties (such as pulsar pe- 
riod and magnetic field derived in the frame of an outer gap 
model), and that the observed y-ray spectrum can be fit to this 
model. Although an unpulsed exc ess from E GRET has been re- 
ported with a significance of 9<x jNel etalJl996h . a significant 
periodicity could not be established. Additionally an ASCA 
X-ray source possibly connected to the EGRET data above 



1 GeV jRoberts. Romani & KawafeOOlh was found in this re- 
gion. Recentlv. lNolanetal] J2003I) reassessed the variability of 
the EGRET source and found a weak variability, which led the 
authors to consider the source finally as a PWN candidate in 
the EGRET high-energy y-ray energy range above 100 MeV. 

Here we report on re-observations of the VHE y-ray source 
HESS J 1825-1 37 and the region surrounding PSR J 1826-1 334 
performed with H.E.S.S. in 2005. H.E.S.S. consists of four 
imaging atmospheric Cherenkov telescopes and detects the 
faint Cherenkov light from y-ray induced air showers in the 
atmosphere above an energy threshold of 100 GeV up to sev- 
eral tens of TeV. Each telescope i s equipped with a mirror area 
of 107 m 2 JBernlohret alJl2003l) and a 960 photo-multiplier 
camera for the detection of the faint Cherenkov light. The tele- 
scopes are operated in a coincidence mode in wh ich at least 
two t elescopes must have triggered in each event jFunk et all 
120041) . The H.E.S.S. system has a point source sensitivity above 
100 GeV of < 2.0 x 10~ 13 cnrV 1 (1% of the flux from the 
Crab nebula) for a 5cr detection in a 25 hour observati on. The 
system is located in the Khomas Highland of Namibia JHintorJ 
2004) and began operation in December 2003. 

2. H.E.S.S. observations of PSR J1 826-1 334 

First indications of a VHE y-ray signal in the region sur- 
rounding the pulsar P SRJ1826-1334 during the H.E .S.S. 
Galactic plane survey JAharonian et alJ l2005bl l2Q06ch trig- 
gered pointed re-observations of the region, resulting in the 
dete ction of an 8.1 <x significa nce signal - named HESS J 1825- 
137 JAharonian et alJl2005cl) . This significance was obtained 
using events within a circle of a priori chosen radius 9 = 0.22° 
from the best fit position as used in a blind search for somewhat 
extended sources. Using a larger integration radius of 9 = 0.4°, 
appropriate to contain most of the emission region, the signif- 
icance increased to 13.4cr. HESS J1825-137 was reported to 
extend -1° asymmetrically to the south of PSR J 1826-1 334 
and shows the same asymmetric exte nsion as the X-ray PWN 
G 18.0-0.7 on a much larger scale. In lAharonian et all J2005cl) 
an association has been proposed between the H.E.S.S. source 
and the X-ray emission region. The different sizes in the two 
energy bands were explained by the difference in the syn- 
chrotron cooling lifetimes of the (higher energy) X-ray emit- 
ting and the (lower energy) IC-y-ray emitting electrons. The 
energy spectrum of the source in the 2004 data within the larger 
integration circle of 9 = 0.4° was fitted by a power law of 
photon index 2.40 + 0.09 sta t ± 0.2 sys at a flux level correspond- 
ing to 20% of the flux from the Crab nebula above 1 TeV. The 
peak of HESS J 1825-1 37 is located just outside the 95% confi- 
dence limits on the position of t he unidentifi ed EGRET source 
3EGJ1826-1302. As shown bv lAharonian et alJ J2005cl) . the 
H.E.S.S. energy spectrum can be connected to the EGRET 
spectrum by extrapolation. Therefore, despite the somewhat 
marginal spatial coincidence, an association between these two 
objects was considered. 

HESSJ1825-137 was revisited in 2005 for ~7 hours in 
pointed observations between June and July and was addi- 
tionally in the field of view of a large part of the pointed 
observations on the nearby (distance ~1°) y-ray emitting mi- 
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Fig. 2. Slices in the uncorrelated excess map of HESS J 1825-1 37 to further illustrate the morphology. The width of the slices is 
0.6° .The direction is given by the orientation from the fit of an elongated Gaussian (see Fig.^ and has a value of 17° (central 
panel) and perpendicular to this direction (right panel). The left panel shows the region in which the slices were taken as black 
boxes overlaid on the H.E.S.S. excess map as shown in Fig.^ The colors denote the 20%, 40%, 60% and 80% intensity contours 
of the VHE y-ray emission. The slices are centred on the best fit position, the position of the pulsar in the slices is marked as 
a dashed black line. As a large part of the observations are taken south of HESS J 1825-1 37, the range beyond 1° north of the 
pulsar is at the edge of the field of view of most observations, and the background subtraction is less reliable in this region. 



croquasar LS5039 (HESS J 1826-14 6), adding another 50. 9 
hours between April and September JAharonian et al.ll2005dh . 
Here we report on the total available dataset (i.e. 2004 and 
2005 data) that includes now ~67 hours of observations with 
HESS J1825-137 within 2.0° of the pointing position of the 
telescopes. The exposure adds up to a total dead-time corrected 
lifetime of 52.1 hours after quality selection of runs accord- 
ing to hardware and weather conditions, thereby increasing the 
observation time by more than a factor of 6 compared to ear- 
lier publications. The mean zenith angle of the dataset pre- 
sented here is 20.1°, the mean offset of the peak position of 
HESS J1825-137 from the pointing direction of the system is 
1.2°. 

The standard H.E.S.S. event reconstruction scheme was ap- 
plied to the raw data a fter calibration and tail -cuts cleaning 
of the camera images d Aharonian et al.ir2004al) . The shower 
geometry was reconstructed based on the intersection of the 
image axes, providing an angular resolution of ~0.1° for in- 
dividual y-rays. Cuts on scaled width and length of the im- 
age (optimised on y-ray simulations and off-source data) are 
applied to s elect y-ray candidates a nd suppress the hadronic 
background JAharonian et aljE004bl) . The energy of the y-ray 
is estimated from the total image intensity taking into account 
the shower geometry. Th e resulting energy resolution is ~ 15%. 
As previously described JAharonian et aljE005bl Eo06dh . two 
sets of quality cuts are applied. For morphological studies of 
a source a rather tight image size cut of 200 photo-electrons 
(p.e.) is applied (along with a slightly tighter cut on the mean 
scaled width), yielding a maximum signal-to-noise ratio for a 
hard-spectrum source. For spectral studies the image size cut is 
loosened to 80 p.e. to extend the energy spectra to lower ener- 
gies. Different metho ds are applied to derive a back ground es- 
timate as described bv lHinton. Berge & FunkU2005i) . For mor- 
phological studies the background at each test position in the 
sky is either derived from a ring surrounding this test position 
(with radius 1 .0°, an area 7 times that of the on-source area, tak- 



ing into account the changing acceptance on the ring), or from 
a model of the system acceptance, derived from off-data (data 
with no y-ray source in the field of view) with similar zenith an- 
gle. In all background methods, known y-ray emitting regions 
are excluded from the background regions to avoid y-ray con- 
tamination of the background estimate. All results presented 
here have been obtained consistently with different background 
estimation techniques. 

3. VHE y-ray emission from HESS J1 825-1 37 

To illustrate the overall morphology of HESS J 1825-1 37, 
Fig- n shows a smoothed excess map of the field of view sur- 
rounding the source, corrected for the changing relative ac- 
ceptance in the field of view. The background for this map 
has been derived from a model of the system acceptance ob- 
ta ined from off-data (sim ilar to the background estimation 
in lAharonian et alJ J2006dl) ). The map has been smoothed with 
a Gaussian of width 2.5'. The inset in the bottom left corner 
shows a Monte-Carlo simulated point-source as it would ap- 
pear in the same dataset taking the smoothing and the point- 
spread function (PSF) for this dataset into account. The pul- 
sar PSR J1826— 1334 is marked by a white triangle. To the 
south of HESS J 1825-1 37, another VHE y-ray source, the 
point -source microquasar LS 5039 (HESS J 1826-148), is vis- 
ible JAharonian et al.l l2005dl> . The color scale for this latter 
source is truncated and thus its apparent size is exaggerated. 
Also shown in Fig. ^ is the 95% positional confidence con- 
tour of the unidentified EGRET source 3EG J1826-1302 (dot- 
ted white), that is possibly associated to HESS J1825-137. 

HESS J1825-137 shows a clearly extended morphology 
with respect to the PSF, extending to the south-west of the pul- 
sar. The position and extension of HESS J 1825-1 37 have been 
determined by fitting the uncorrelated (i.e. unsmoothed) excess 
map to a model of a 2-D Gaussian y-ray brightness profile of 
the form p oc exp(-0 2 /2cr; ource ), convolved with the PSF for 
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Fig. 3. Energy flux E 2 dN/dE of HESS J1825-137 in y-rays as measured by H.E.S.S. (black filled circles) up to -35 TeV for 
a large part of the emission region within an integration radius of 0.8° around the best fit position. Also shown as open circles 
are the energy flux points of the possibly related unidentified EGRET source 3 EG J 1826-1 302. The background estimate for the 
spectrum has been derived from off-runs (i.e. runs without y-rays in the field of view). The spectrum can be fit by a power law 
dN/dE = I a E- r but the poor X 2 of the fit (40.4 for 15 de grees of freedom) suggests a different spectral shape. Here a power law 
fit with an exponential cutoff is shown, but it should be noted that this shape is not favoured over other similarly complex spectral 
shapes as discussed in the text in more detail. 



this dataset (68% containment radius: 0.075°). The best fit po- 
sition - equivalent to the center of gravity of the source - is at 
18h25m41s +3 stat s , -13°50'21"+35 stat "(here and in the fol- 
lowing the epoch J2000 is used), the best fit rms extension 
is o- source = 0.24° + 0.02 stat o . However, the x 2 P er degree 
of freedom is not satisfactory (1295/1085), indicative of the 
more complex morphology of the source. Reflecting the non- 
Gaussian and skewed source profile, the position of the peak 
in the y-ray emission (at 18h25m57s, -13°43'36.8"as deter- 
mined by fitting a 2-D Gaussian in a restricted region around 
the peak) is slightly shifted at a distance of ~8' to the best fit 
position. The pulsar PSR J 1826-1 334 is located at a distance 
of ~10' from the peak y-ray emission and ~17' from the best 
fit position. To test for a different source morphology, an elon- 
gated Gaussian with independent cr source along the minor and 
^source a l° n g the major axis and a free position angle u> (mea- 
sured counter-clockwise from the North) has also been fitted. 
This elongated fit gives a best fit position of 18h25m41s ±4 stat s, 
-13°50'20"+40 sta t", consistent within errors to the symmetri- 
cal fit. The fit yields only a slight indication for an elongation 
with cr source = 0.23 +0.02 stat and a-' som . ce = 0.26 +0.02 stat at a po- 
sition angle of a> — 17° + 12° stat . The;^ 2 per degree of freedom 
(1288/1083) is still relatively poor. The best fit position devi- 
ates slightly from the best fit po sition reported in earlier papers 
( I Aharonian et al.l2005blE006cl) . The difference can mainly be 
attributed to the different fit range. The best fit parameters of the 
elliptical fit are shown as a black square and ellipse in Fig.^ 
Note that the fitted position angle is consistent within errors 



with the orientation of the line connecting the pulsar position 
and the best fit position, which amounts to 23.1°. 

Fig. |2 shows slices in the direction of the position angle 
(17°) of the elliptical fit (centre) and in the direction perpendic- 
ular to it (right). The width of the slices is chosen to be 0.6°, 
the slices are illustrated in the left panel as black dashed boxes. 
The position of the pulsar in the slices is marked as a dashed 
black line. It can be seen, that the peak of the H.E.S.S. emis- 
sion is close to the pulsar position but slightly shifted as is also 
apparent from the two-dimensional excess plot. Also visible in 
the central panel is the rather sharp drop from the peak posi- 
tion towards the north-eastern direction and the longer tail to 
the south-western direction. Some indication for an additional 
excess to the north of HESS J1825-137 is seen in Fig.^and 
in the central panel of Fig. 13 at a distance of ~0.7° from the 
pulsar position. Further investigation of this feature will have 
to await future data, in particular given that most current data 
were taken on positions south of the pulsar, with regions in the 
north near the edge of the field of view. 

For the spectral analysis the image size cut is loosened to 
80 p.e. to achieve a maximum coverage in energy. The result- 
ing spectral analysis threshold for the dataset described here is 
270 GeV. Events with reconstructed direction within an angle 
6 = 0.8° of the source location are considered on-source. No 
correction for the y-ray emission extending beyond this angu- 
lar cut has been applied. Thus the flux level determined cor- 
responds to the flux level of the source within the integration 
region and might be an underestimation of the flux from the 
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Table 1. Fit results for different spectral models for the whole emission region within an integration radius of 0.8° around 
the best fit position and the background derived from off-data. The differential flux normalisation /q is given in units of 
10~ 12 cirT 2 s" 1 TeV~'. E, E B , and E c are given in units of TeV. The last column gives the integrated flux above the spectral 
analysis threshold of 270 GeV in units of 10~ n cirT 2 s _1 . The power-law fit provides a rather poor description of the data. Thus 
fits of a power law with an exponential cutoff (row 2), a power law with an energy dependent photon index (row 3), and a broken 
power law (row 4; in the formula, the parameter S — 0.1 describes the sharpness of the transition from Ti to T2 and is fixed in 
the fit) are also given. Note that some of the fit parameters are highly correlated. 



Fit Formula for ^ 


Fit Parameters 


X 2 (ndf) 


F1UX >2 7() GeV 


kE- 1 

Iq E~ t exp(-E/E c ) 

k E- r+l1 1oe£ 

/„ {E/E B Y r > (l+(E/E B ) lls ) Sir, ~ F2) 


7 = 19.8 ±0.4 r = 2.38 + 0.02 

I a = 21.0 ±0.5 T = 2.26 + 0.03 E c = 24.8 ± 7.2 

I a = 21.0 ±0.4 T = 2.29 ±0.02 f} = -0.17 ± 0.04 

/o = 2.2 ±1.0 r,= 2.26 ±0.03 T 2 = 2.63 ± 0.07 E B = 2.7 ± 0.5 


40.4 (15) 
16.9 (14) 

14.5 (14) 

15.1 (13) 


87.4 ± 2.0 

86.7 ± 2.5 

82.8 ± 2.2 
84.6 ± 38.5 



whole source. In the determination of the energy spectrum, the 
energy of each event is corrected for the time-varying detector 
optical efficiency, relative to that used in Monte Carlo simula- 
tions to estimate the effective area of the instrument. The opti- 
cal efficiency is estimated fro m single muon events detecte d 
during each observation run llLerov et alJl2003t lBoldl20ol . 
The mean energy correction is ~25%. For the spectral analy- 
sis the background is taken from positions in the same field of 
view with the same offset from the pointing direction as the 
source region. This approach is taken to avoid systematic ef- 
fects from the energy-dependent system acceptance function 
(which is to a good approximation radially symmetric). In an- 
other approach off-data have been used in the background esti- 
mation to confirm the results from the same field of view, using 
either the same shaped region as the on-region in the off-data 
or using again off-regions distributed with the same offset from 
the pointing direction of the system as the on-region. The total 
significance of the emission region with the loose cuts is 33.8<r 
with an excess of 19510 + 577 y-ray events. Fig.[3]shows the 
spectral energy distribution in terms of energy flux E 2 dN/dE 
of the H.E.S.S. emission region (full black circles). Also shown 
are the energy flux points and the spectral fits of the possibly 
related unidentified EGRET source 3 EG J 1826-1 302 (open cir- 
cles). Given the poor angular resolution of EGRET, these data 
are taken on a scale similar to that of the full H.E.S.S. emis- 
sion region and can thus be compared to the total H.E.S.S. 
flux. From this figure one can see that the unidentified EGRET 
source 3EG J1826— 1302 could be associated with the H.E.S.S. 
emission region from a spectral continuity point of view. 

A fit of the differential energy spectrum from 270 GeV 
up to ~35 TeV by a power law dN/dE = Iq E~ r yields a 
normalisation of Iq = 19.8 + 0.4 stat + 4.0 sys x 10~ 12 TeV" 1 
cm" 2 s _1 and a photon index F = 2.38 + 0.02 stat + 0.15 sys 
(see Tabled The flux of HESS J1825-137 above 1 TeV corre- 
sponds to ~68% of the flux from the Crab nebula. Note that 
this flux is significantly hig her than the previously reported 
flux ( I Aharonian et al.ll2005d) due to a significantly increased 
integration radius (0.8° instead of 0.4°) in the attempt to cover 
the whole source region. Integrating only within the smaller 
region of 0.4° the flux level is consistent with the previously 
published result. The power-law fit represents a rather bad de- 
scription of the data (as can be seen x 1 of the fit) and suggests 
therefore a different spectral shape. Various models have been 



fit to the data to investigate the shape of the spectrum. Tabled 
summarises these fits. Three alternative shapes have been used: 
a power law with an exponential cutoff E c (row 2), a power law 
with an energy dependent exponent (row 3), and and a broken 
power law (row 4). In all cases, Iq is the differential flux nor- 
malisation, and the photon indices are specified as F. It is ev- 
ident that the alternative descriptions of the spectrum describe 
the data significantly better than the pure power law as can be 
seen from the decreasing^ 2 /ndf (see TableQ. 

Given the large dataset with more than 19,000 y-ray ex- 
cess events and given the extension of HESS J 1825-1 37, a spa- 
tially resolved spectral analysis has been performed to search 
for a change in photon index across the source, similar to the 
detailed a nalysis of the y-ray SNR RX J 17 13.7-3946 as per- 
formed in lAharonian et al] J2006dl) . Fig.|4]shows energy spec- 
tra determined in radial bins around the pulsar position, cov- 
ering the extended tail of the VHE y-ray source. The inset of 
Fig.|4]shows again the H.E.S.S. excess map as shown in Fig.^ 
along with wedges that illustrate the regions in which the en- 
ergy spectra were determined, with radii increasing in steps of 
0.1°; the innermost region is centred on the pulsar PSR J 1826- 
1334. The opening angle of the wedges was constrained by 
LS 5039 in the southern part and by the apparent end of the 
emission region in the northern part. For all regions the energy 
spectrum has been determined by defining the wedge as the on- 
region. The background estimate has been derived from circles 
distributed on a ring around the pointing direction. The radius 
of this ring was chosen to be equal to the distance of the centre 
of gravity of the wedge to the pointing direction. This approach 
ensures a similar offset distribution in the on- and off-dataset 
and has been used to determine the background estimate from 
the same field of view as well as from off-data taken on regions 
without y-ray sources. Consistent results were achieved in both 
methods. 

Along with each spectrum in Fig. |4] the power law fit to 
the innermost region centred on the pulsar position is shown as 
a dashed line for comparison. A softening of the energy spec- 
tra is apparent with increasing distance from the pulsar. This 
softening is equivalent to a decrease of the source size with in- 
creasing energy and provides the first evidence for an energy 
dependent morphology detected in VHE y-rays. Differences in 
the energy bin sizes arise from the fact that for non-significant 
photon points the bin size was increased. It has been verified 
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Fig. 4. Energy spectra in radial bins. Inset: H.E.S.S. excess map as shown in Fig.^ The wedges show the radial regions with 
radii in steps of 0.1° in which the energy spectra were determined. The innermost region is centred on the pulsar PSR J 1826- 
1334. Main Figure: Differential energy spectra for the regions illustrated in the inset, scaled by powers of 10 for the purpose of 
viewing. The spectrum for the analysis at the pulsar position is shown as a reference along with the other spectra as dashed line. 
For all regions the energy spectrum has been determined as described in the text and has been fitted by a power-law in a restricted 
energy range between 0.25 and 10 TeV. 



that this approach does not change the result of the fit. Due 
to the different distribution of offsets from the pointing direc- 
tion of the system in the different regions, the photon analysis 
threshold changes slightly, thus some of the different spectra do 
not start at exactly the same energy. 

Fig-IUsummarises the findings of Fig.|4]by plotting the fit 
parameters of the power law fit versus the distance of the region 
to the pulsar position. Shown are the results using two differ- 
ent background estimation techniques in the spectral analysis. 
The left panel shows the photon index as a function of the dis- 
tance from the pulsar. A clear increase of the photon index for 
larger distances from the pulsar position is apparent; the pho- 
ton index seems to level off within errors to a value of ~2A 
+0.1 at a distance of ~0.6°. The right panel shows the surface 



brightness (i.e. the integrated energy flux EdN/dE per unit area 
between 0.25 TeV and 10 TeV) as a function of the distance to 
the pulsar position. Again here it can be seen, that the maxi- 
mum of the emission is slightly shifted away from the pulsar 
position as was already apparent in Fig.|2] In both panels, the 
error bars denote ±lcr statistical errors. Systematic errors of 
20% on the flux and 0.15 on the photon index are to be as- 
signed to each data point in addition. However, since all spectra 
come from the same set of observations, these systematic errors 
should be strongly correlated, and will cancel to a large extent 
when different wedges are compared. Table |2] summarises the 
different spectral parameters determined in the wedges using 
the reflected background from the same field of view. 
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Fig. 5. Energy spectra in radial bins. Left: Power-law photon index as a function of the radius of the region (with respect to the 
pulsar position) for the regions given in Fig.|4] The closed points are obtained by deriving the background estimate from regions 
with the same offset as the on-region within the same field of view. The open points are derived using off-data (data without y-ray 
sources in the field of view) for the background estimate. A constant fit to the filled points yields ax 2 of 58 for 11 degrees of 
freedom, clearly showing the dependence of the photon index on the distance from the pulsar position. Right: Surface brightness 
between 0.25 and 10 TeV per integration region area in units of 10~ 12 erg crrT 2 s" 1 deg~ 2 as a function of the distance to the pulsar 
position. Error bars denote + \cr statistical errors. 



Table 2. Spectral parameters for the radial bins surrounding PSR J1826-1334. PSF denotes a H.E.S.S. point-source analysis at 
the pulsar position. The background estimate for the numbers in the table have been derived from reflected positions within the 
same field of view. The energy flux and surface brightness are given for the energy range between 0.25 and 10 TeV. 



Radius 


Photon Index 


*7ndf 


Area 


Energy Flux 


Surface Brightness 


(deg) 








(deg 2 ) 


(erg cm 


2 S" 1 ) 


(erg cirT 2 s deg -2 ) 


PSF 


1.83 


±0.09 


2.8/3 


3.1e-02 


2.9e-12 ± 


3.2e-13 


9.3e-ll ± 1.0e-ll 


0.15 


1.96 


+ 0.08 


0.8/3 


2.4e-02 


2.9e-12 ± 


2.4e-13 


1.2e-10± 1.0e-ll 


0.25 


2.20 


±0.06 


3.1/3 


3.9e-02 


4.3e-12 ± 


2.5e-13 


l.le-10±6.4e-12 


0.35 


2.25 


±0.06 


6.9/3 


5.5e-02 


5.3e-12 ± 


2.8e-13 


9.7e-ll ± 5.1e-12 


0.45 


2.32 


±0.07 


7.1/3 


7.1e-02 


6.2e-12 ± 


3.7e-13 


8.8e-ll ± 5.2e-12 


0.55 


2.37 


±0.06 


8.5/3 


8.6e-02 


6.9e-12 ± 


3.9e-13 


7.9e-ll ± 4.5e-12 


0.65 


2.36 


±0.08 


0.4/3 


1.0e-01 


5.7e-12 ± 


3.6e-13 


5.5e-ll ± 3.5e-12 


0.75 


2.41 


±0.09 


8.3/3 


1.2e-01 


5.4e-12 ± 


4.0e-13 


4.6e-ll ± 3.4e-12 


0.85 


2.42 


±0.09 


6.0/3 


1.3e-01 


5.1e-12 ± 


4.2e-13 


3.8e-ll ± 3.1e-12 


0.95 


2.59 


±0.13 


2.4/3 


1.5e-01 


5.2e-12 ± 


5.4e-13 


3.5e-ll ± 3.6e-12 


1.05 


2.43 


±0.09 


6.4/3 


1.6e-01 


5.9e-12 ± 


4.9e-13 


3.6e-ll ± 3.0e-12 


1.15 


2.45 


±0.35 


3.4/3 


1.8e-01 


2.9e-12 ± 


6.8e-13 


1.6e-ll ± 3.8e-12 



Whereas the H.E.S.S. observation of an energy depen- 
dent morphology represents the first detection of such an ef- 
fect in y-ray astronomy, the dependence of observed pho- 
ton index on radius (commonly known as "T - r" relation) 
is well known from X-ra y observations of P WN other than 
the Crab. For G2 1.5-0.9 Isiane et alJ d2000h found a value 
of F ~ 1.5 near the PWN termination shock, after which 
it converges to a value of ~2.2 in the outer nebula. For 
the PWN 3C58, the photon index in creases from 1.9 i n the 
torus to - 2.5 at the edge of 3C 58 dSlane et al.ll2004 . For 
this object iBocchino et alJ ll200ll) . also found that the energy 
flux per radial interval for 3C58 remains approximately con- 
stant, consistent with our findings for HESS J 1825-1 37. A 
similar constant energy flux with increasing radius was also 



found in X-ray observations of another VHE r-r ay source 
GO. 9+0.1 JPorquet. Decourchelle & Warwicld Eooi." For this 
composite remnant, the photon index also varies with ra- 
dius from 1.5 (beyond the compact core) to ~2.5 near the 
edge of the PWN. In the case of the more evolved Vela 
PWN iMangano et alJ J2005I) found a radial variation of 1.55 
to 2.0. For most of these remnants a total change in the photon 
index of ~0.5 is seen, as expected for cooling losses. Attempts 
to model the T — r relationship wer e not successful in the past - 
ISlane et alJ d2004 showed that the lKennel & Coronitil ( ll984bl) 
model for convective flow (which includes the conservation of 
magnetic flux) fails to reproduce this well-known T — r rela- 
tionship for PWN which are evolved beyond the Crab phase. It 
should additionally be noted, that the H.E.S.S. observation of 
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an energy dependent morphology is the first unambiguous de- 
tection of a spectral steepening away from the pulsar, for fixed 
electron energies; in X-rays the situation is complicated by a 
possible variation of the magnetic field with increasing distance 
from the pulsar; if the X-ray spectrum is probed near or above 
the peak of the SED, a variation of the field will influence the 
slope. Depending on the age and magnetic field, one might ex- 
pects to see similar effects in other VHE y-ray PWN, but so far 
only HESS J 1825-1 37 has sufficient statistics to clearly reveal 
the energy dependent morphology. 

To further investigate the spectral properties of 
HESS J1825-137, the emission region has been segmented 
into boxes. The result of the spectral analysis in these boxes is 
shown in Fig. [6] The left panel shows in red VHE y-ray excess 
contours as given in Fig.^ Overlaid are 12 boxes for which 
spectra were obtained independently. The photon index result- 
ing from a power law fit in each region is grey-scale coded in 
bins of 0. 1 . Also here a softening of the spectral indices away 
from the pulsar position is apparent, although the error bars are 
larger than in Fig.[5]due to the smaller integration regions. The 
size of the boxes is equivalent to the ones use d in the analysis 
of the shell-type SNR RX J1713.7-3946 JAharonian et all 
l2006dl) . where no spectral variation has been detected. The 
right hand figure shows the correlation of photon index T to 
integral flux per square degree above 1 TeV. A mild correlation 
between the flux per deg 2 and the spectral index exists and 
the correlation coefficient between these two quantities is 
-0.46 + 0.14. 

4. Interpretation 

Obviously, the key aspect of the new H.E.S.S. data is the soft- 
ening of the energy spectrum at larger distances from the pul- 
sar position or equivalently, the decrease in source size with 
increasing energy of the y-rays. The following discussion is 
based on the plausible assumption that particles have been 
accelerated in the proximity of the pulsar during the last 20 
kyears - the characteristic age of the pulsar - and that they 
were then transported by diffusion and/or convection to form 
the extended nebul a of about 10-15 pc ra dius (for a recent 
review of PWN see lGaensler & Slanell2006|) . The asymmetric 
shape of the nebula most likely reflects the ambient conditions, 
and could be caused by a reverse supe rnova shock created in 
the dense medium north of th e pulsar Jvan der Swaluw & Wul 
1200 lU Aharonian et alJl2005ct) . geometrically by the motion of 
the pulsar leaving a PWN 'trail', or even dynamically by a bow 
shock, resulting in a 'cometary' tail; however, the high speeds 
in excess of 1000 km/s required in case of HESS J1825-137 
favour the first scenario. Spectral variation with distance from 
the pulsar could result from (i) energy loss of particles dur- 
ing propagation, with radiative cooling of electrons propagat- 
ing outward from the pulsar termination shock as the obvious 
loss mechanism, as well as from (ii) energy dependent diffu- 
sion or convection speeds, and from (iii) variation of the shape 
of the injection spectrum with age of the pulsar which, after 
propagation, translates into a spatial variation of spectra. 

Loss mechanisms in (i) include, e.g, adiabatic expansion, 
ionisation loss, bremsstrahlung, synchrotron losses and inverse 



Compton (IC) losses; only the last two result in a lifetime 
r = E/(dE/dt) which decreases with energy and hence causes 
power-law spectra to steepen, due to the quadratic dependence 
of dE/dt on the particle energy jBlumenthal & Gouldll970l 
lKardashevl[l96l . A source size which decreases with energy 
is therefore a strong indication that the accelerated particles are 
electrons. The lifetime due to synchrotron and IC losses is: 

Here, w r = J]w p h + 0.025 eVcirT 3 (B//iG) 2 is the total en- 
ergy density in the form of radiation and magnetic fields, 77 
is a normalisation factor to account for the reduction of IC 
losses due to the transition into the Klein-Nishina regime and 
B the magnetic field. Given the density of the 2.7 K CMBR 
W2.7K = 0.26 eVcirT 3 , and the fact that IC losses of multi- 
TeV electrons on the diffuse optical/IR photons are strongly 
suppressed due to the Klein-Nishina effect, and even for the 
CMBR are reduced by a factor rj » 2/3 at H.E.S.S. energies, 
synchrotron losses dominate for B > 3 /uG. In case of continu- 
ous injection and radiative lifetimes short compared to the age 
of the source, T ia d(E) <& T, the spectral index a of the elec- 
trons steepens by one unit, corresponding to a change of the 
photon index by half a unit, which approximately matches the 
observed variation between the inner and outer regions of the 
nebula (Fig. O- In the Thomson regime, the energy of the par- 
ent electrons is E e « 20(£y/TeV ) 1/2 TeV, corresponding to 
the range from about 10 to 100 TeV for y-ray energies between 
0.2 and 20 TeV. Cooling time scales below 20 kyears require 
B > 6 juG for 10 TeV electrons; at 100 TeV the lifetime is 
below 20 kyears already in typical 3 juG interstellar fields, so 
some steepening of spectra at the highest energies is expected 
even in relatively modest fields. 

It is then instructive to consider the energy budget of the 
PWN in an electronic scenario. The assumed large distance 
of a* 4 kpc and the relatively high y-ray flux, F - 1.5 x 
10~ 10 erg/cm 2 s above 200 GeV, imply a quite luminous VHE 
y-ray source, L y ~ 3 x 10 35 erg/s. This luminosity is compara- 
ble to that of the Crab nebula, while the spin-down luminosity 
of the pulsar is smaller by two orders of magnitude. Thus, the 
efficiency of the y-ray production in HESS J 1825-1 37 is much 
higher, e y = L y (> 200Ge V)/L rnt =a 0.1. A relatively large effi- 
ciency is not unexpected jAharoniaa Atovan & Kifundl 19971) 
since the much lower magnetic field in a nebula powered by a 
less energetic pulsar results in a more favourable sharing be- 
tween IC and synchrotron energy losses. In a steady state, and 
neglecting non-radiative energy losses, the efficiency for y-ray 
production is 

e Trad ~. ^ (2) 

e y ~ e e ^ « e e ^ + Q 025eVcm -3(B/ / uG) 2 

where tic is the lifetime due to IC losses and e e is the fraction 
of pulsar spin-down power going into 10 to 100 TeV electrons, 
corresponding to the observed y-ray energies. Unless the elec- 
tron spectral index is well below 2 at energies below 10 TeV, e e 
will be below 10-15%, taking into account the sharing of spin- 
down energy between particle and field energies. An efficiency 
e y of 0. 1 cannot be obtained, even for rather small magnetic 



10 



F. Aharonian et al.: Energy dependent y-ray morphology in HESS J1825— 137 




18 h 30 m 18 h 20 m 2 2.2 2.4 2.6 2.8 

RA J2000 (hours) Photon index r 



Fig. 6. Spatially resolved spectral analysis of HESS J1825— 137. Left: Shown in solid red are VHE y-ray excess contours from 
Fig-ID linearly spaced at the 20%, 40%, 60% and 80% maximum intensity levels. Superimposed are the 12 boxes (0.26° x 0.26°) 
for which spectra were obtained. The black dashed line is the 0.8° radius circle that was used to integrate events to produce a 
spectrum of the whole region. The photon index obtained from a power-law fit in each region is colour coded in bins of 0.1. 
Right: Integral flux above 1 TeV per square degree versus photon index, for the 12 regions shown in the left panel. The error bars 
are ±l<x statistical errors. The correlation coefficient for the two quantities is -0.46 + 0.14. 



fields in the range of a few /uG. Detailed numerical simula- 
tions with (optionally time-dependent) electron injection and 
cooling confirm that an energy input about one order of magni- 
tude higher than the current spin-down luminosity is required 
to sustain the observed gamma-ray flux and to quantitatively re- 
produce the measured spectrum, assuming that the distance of 
~ 4 kpc is correct. A likely solution is that the spin-down power 
of the pulsar was significantly higher in the past; for modest B 
fields of a few fiG electron lifetimes in particular at lower ener- 
gies are of the order of the pulsar age and the time variation of 
spin-down luminosity needs to be taken into account. For ex- 
ample, with L lot oc r 2 for a braking index of n — 3, 'relic' elec- 
trons released in the early history of the pulsar and surviving 
until today can provide sufficient energy. To allow accumula- 
tion of electrons over the history of the pulsar, magnetic fields 
should not exceed 10yuG. 

A discussion of the energy-dependent morphology requires 
assumptions concerning the transport mechanism. At least in 
the inner regions of the nebula, convection is likely to domi- 
nate over diffusion. Indeed, the variation of surface brightness 
across the source - roughly proportional to 1/9, where 9 is the 
angular distance from the pulsar (see Table f2J - is difficult to 
account for in purely diffusive propagation. A surface bright- 
ness oc Q~" is obtained - for spherical symmetry - from a vol- 
ume density oc r~ n . Neglecting cooling effects, a 1/9 depen- 
dence is hence obtained for a constant radial convection veloc- 
ity, resulting in a 1 jr 2 density distribution. For constant con- 
vection speed, energy conservation requires a rapid decrease of 
fi-fields with distance from the pulsar, with very low fields at 
the edge of the PWN unless one is dealing with a very strong 
and yo ung source such as the Crab nebula ( iKennel & Coroniti 
Il984bl) . A convection speed v(r) oc 1/r would allow a con- 
stant B-field. Such convection results in constant surface den- 
sity; however, the electron density at a fixed electron energy 



- and therefore the y-ray intensity - will again decrease with 
distance once cooling is included. A speed v(r) oc 1/r results 
in a propagation time t ~ r 2 and, at energies where the elec- 
tron lifetime r ra d oc l/E e is shorter than the lifetime T of the 
accelerator, in a source size R oc E? 1 ^ 2 . 

A similar result is obtained for the diffusion case (ii), which 
is expected to be relevant near the outer edge of the nebula. The 
diffusive source size is governed by the diffusion coefficient 
D(E), which is frequently parametrised in a power-law form 
D(E) = D(){E I Eq) s , with 5 between for energy-independent 
diffusion and 1 for Bohm diffusion. The resulting size can 
be estimated to R [2D(E)t] 1 ^ 2 with the propagation time 
t again given by the age T of the accelerator or the lifetime 
^ad of radiating particles, whatever is smaller. For lifetimes 
^"rad <K T short compared to the age of the accelerator, one ob- 
tains R oc Ef~ . In case of Bohm-type diffusion with 6=1, 
the radiative losses and the diffusion effects compensate each 
other and the size becomes effectively energy independent. For 
energy independent diffusion, i.e. 6 = 0, the size decreases with 
energy again as R oc E"^ 1 . 

Option (iii) - a time-variable acceleration spectrum - is 
a distinct possibility in particular for accelerated electrons. 
Higher pulsar spin-down luminosity in the past will have been 
associated with higher B fields and a lower cutoff energy, gov- 
erned by the relation between acceleration and radiative cool- 
ing time scales. In either case (i), (ii) or (iii), the new H.E.S.S. 
results therefore provide evidence of an electronic origin of 
the VHE y-ray emission, and require that characteristic cool- 
ing time scales are, or at some earlier time were, shorter than 
the age of the nebula. 
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4.1. Conclusion 

We have presented detailed morphological and spectral studies 
of the VHE y-ray source HESS J 1825-1 37 that has been origi- 
nally detected in the survey of the inner Galaxy, conducted by 
H.E.S.S. in 2004. The y-ray spectrum of the source has been 
measured over more than two decades between ~270 GeV and 
~35 TeV. The energy spectrum shows indications for a devia- 
tion from a pure power-law. Several spectral shapes have been 
applied to fit the data and it seems, that a broken power-law 
or a power-law with energy dependent photon index provide a 
better description than a pure power law. The large data set 
has provided the possibility for a spatially resolved spectral 
study. A significant softening of the y-ray spectrum away from 
the position of the energetic pulsar PSR J1826— 1334 has been 
found, providing the first direct evidence of an energy depen- 
dent morphology in VHE y-rays. The studies performed here 
significantly strengthen the case that the VHE y-ray emission 
originates in the wind nebula of PSR J1826-1334. It is diffi- 
cult to explain the measured y-ray luminosity in terms of the 
current spin-down luminosity of the pulsar. A like scenario is 
a significant contribution of 'relic' electrons released in the 
early history of the pulsar, when the spin-down luminosity is 
higher. The variation of index with distance from the pulsar is 
attributed both to IC and synchrotron cooling of the continu- 
ously accelerated electrons. 
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